Direct estimation of decoherence rates 
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The decoherence rate is a nonlinear channel parameter that describes quantitatively the decay of 
the off-diagonal elements of a density operator in the decoherence basis. We address the question 
of how to experimentally access such a nonlinear parameter directly without the need of complete 
process tomography. In particular, we design a simple experiment working with two copies of the 
channel, in which the registered mean value of a two-valued measurement directly determines the 
value of the average decoherence rate. No prior knowledge of the decoherence basis is required. 



I. INTRODUCTION 

Complete quantum tomography is a very complex task 
which rapidly becomes computationally intractable as 
the number of subsystems composing the system un- 
der investigation increases (see e.g., Ref. [1]). More 
specifically, if we assume that d is the dimension of each 
of the subsystems, then a system composed of N such 
subsystems has dimension D — d N . The goal of state 
(or channel) tomography is then to determine D 2 — 1 
or D 2 (D 2 — 1) real parameters, respectively. That is, 
the number of parameters grows exponentially with the 
number N of quantum subsystems that comprise the 
complete system. As a consequence, already for few- 
qubit systems the tomography of associated quantum 
devices requires an almost insurmountable experimental 
and computational effort. 

This enormous effort is the reason why physicists put 
forth tasks that are more modest than complete tomog- 
raphy. For instance, quite often the system under con- 
sideration is not completely unknown; rather, some vital 
prior information is available. If so, the estimation prob- 
lem can become tractable. This idea has led to the de- 
velopment of biased estimation schemes [2-5] that work 
efficiently (i.e., with resource requirements scaling only 
polynomially with system size) for states within a cer- 
tain class of "expected states". Other research lines do 
without estimating all system parameters and properties 
and instead focus on identifying parameters that are ex- 
perimentally feasible, i.e., quantities that can be "easily" 
measured [6, 7]. Here we focus on one such parameter, 
the (average) decoherence rate, and discuss the minimal 
(informationally incomplete) resources needed for its es- 
timation. 

The phenomenon of decoherence is often recognized as 
the main obstacle to the experimental implementation of 
quantum information technologies. In its essence, it con- 
verts a quantum superposition (with associated ampli- 
tudes) into an incoherent mixture (with associated prob- 
abilities). Mathematically, a loss of quantum coherence 
is reflected by a decrease of the off-diagonal terms (ir- 
respective of their initialization) of the system's density 



operator in the so-called decoherence basis, while the di- 
agonal elements remain preserved. In the present paper 
we will assume that the estimated processes are of this 
form, and we will focus on possible methods of estima- 
tion of their parameters, especially of the so-called de- 
coherence rates, defined as fractions of absolute values 
of the final and initial values of off-diagonal elements. 
These parameters are not linear, which opens interesting 
questions on their direct experimental accessibility. Al- 
though in this paper we will consider only a very special 
case of nonlinear parameters, our discussion and findings 
are aiming to develop a general mathematical framework 
dealing with such type of estimation problems. 

In the next section we will briefly describe the relevant 
properties of decoherence channels. In Section III we 
will focus on an experimental setup to access directly 
the decoherence rate (being a nonlinear parameter) for 
qubit channels. We will generalize this consideration to 
the case of d-dimensional quantum systems in Section IV. 
The achieved results will be applied in Section V to direct 
estimation of the (inverse) decoherence rate of a double- 
commutator master equation. In the final Section VI we 
will briefly summarize our results and discuss a possible 
experimental realization. 



II. DECOHERENCE CHANNELS 

Let us denote the Hilbert space of the considered sys- 
tem by W. By £(%) we denote the set of linear operators 
on H, with S(H) being the subset of density operators 
(Trp = 1, q > 0). States are identified with elements 
of <S("H), and (deterministic) quantum processes are de- 
scribed by quantum channels, i.e., by completely positive 
trace-preserving linear maps £ : £(H) —> £(%). These 
properties guarantee that quantum channels map quan- 
tum states into quantum states, i.e., £(S(H)) C S(H). 
A distinguished role within the set of quantum channels 
is played by unitary channels for which £ (g) — UgW, 
where UW = U'U = I. It is quite common to identify 
decoherence channels with non-unitary ones, although for 
most of the channels there is no orthonormal basis that 
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is invariant under the channel's action. If such a basis 
exists we say that the (non-unitary) channel is a pure 
decoherence channel. 

Elementary properties of pure decoherence channels 
were studied in Refs. [8, 9]. Pure decoherence channels 
are parametrized by a choice of the decoherence basis 
and a choice of (complex) inverse decoherence rates, the 
latter independently for each (mutually conjugated) pair 
of off-diagonal terms. In particular, under the action of 
a pure decoherence channel, elements of the density op- 
erator expressed in the decoherence basis undergo the 
transformation 

Qjk Qj k = UjkQjk , (2.1) 

where \oJjk\ are the inverse decoherence rates (j ^ k) 
and, in general, ujjk = uJkj, l w jfe| < 1 and ^jj = 1; 
thus the transformation is characterized by D(D — 1) 
real parameters. Moreover, we need to specify addi- 
tional D(D — 1) parameters in order to determine the 
(unordered) decoherence basis. Altogether, we see that a 
pure decoherence channel is characterized by 2D(D — 1) 
real parameters, which is significantly less than for a 
general channel but still exponential in the number of 
the system's constituents, even if the decoherence ba- 
sis is known. Our goal is to design an experiment that 
would allow us to determine the value of just a single pa- 
rameter which characterizes the "average" decoherence. 
Namely, we will investigate measurement of the quan- 
tity A 2 := rf ( d 2 _ 1 ) X)j<fc l^jfcl 2 , which for qubits reduces 
to estimation of the inverse decoherence rate A. Simulta- 
neously, we require that the experiment provides as little 
information as possible about any other parameters. In 
an ideal case the inverse decoherence rate A will be mea- 
sured as the mean value (or its function) of a specific 
two- valued measurement. If this happens to be the case, 
we say the parameter is measured directly. 



III. QUBIT CASE 

The most general qubit pure decoherence channel can 
be written in the following form [9] 

q^q' = £{q) =pg+(l- P )UgU\ (3.1) 

where U is an arbitrary unitary operator. The eigen- 
basis of U determines the decoherence basis. Suppose 
U = e io |0)(0| + e i6 |l)(l|. Then oj 01 = p + (1 - p)e^ a - b \ 
and the (unique) qubit inverse decoherence rate reads 
A = |cjoi| = ^/cjoi^io- Our goal is to propose an "opti- 
mal" measurement that would allow us to determine this 
parameter directly, while measuring as little redundant 
other information as possible. In order to achieve our 
goal we have to specify an initial probe state u> and a 
two-valued observable (outcomes associated with effects 
E,I — E) such that 

\=f (Tr(£ x ® T){u]E) , (3.2) 



where / is a function on the interval [0,1]. According 
to the Choi-Jamiolkowski isomorphism, linear maps T 
from C(H) to £(/C) are in one-to-one correspondence 
with linear operators F defined on H ® JC. In particular, 

f = (,f®x)[q + ] with n+ = Ej,fc=i \j®j)(k®k\ e n®n 

being an unnormalized version of the corresponding max- 
imally entangled state. The complete positivity of chan- 
nels translates into positivity of the operators. Writing 
co = (X®J-)[0+], il x = (£ A ®X)[ft+] undF = {I®T*)[E\, 
Eq. (3.2) can be rewritten in the form 

A = / (TtVI\F) , (3.3) 

where F is an element of a so-called process POVM 
(PPOVM) (introduced in Refs. [13, 14]) which fully cap- 
tures the adjustable degrees of freedom (choice of initial 
probe state and final measurement) in the experiment. 

A. Known decoherence basis 

Let us suppose that the decoherence basis is known to 
be the one in which the operator Sl + is defined. Then its 
Choi-Jamiolkowski representation reads 

n x = |oo)(oo| + |ii)(ii| + woi|oo)<ii| +wi |ii><ocM) 

It follows from Eq. (3.3) that no single Hcrmitian oper- 
ator Q allows for a determination of A = ^/cjqiWio = 
Tr[f2A<5], simply because this parameter is not linear. 
However, there are several ways to determine the value 
of A by combining expectation values of more than one 
Hermitian operator. Since A = |(ll|fi A |00)| =: |Trfi A T|, 
the decoherence parameter A can be accessed "directly" 
by measuring the mean value of a non-Hermitian opera- 
tor T = 1 00) (11 1. In practice, this requires to realize an 
experiment according to the following recipe: 

1) Prepare a maximally entangled bipartite state and 
send one of its subsystems through the decoherence chan- 
nel. 

2) Define observables X := {A+,A_,A Q } and Y := 
{Y+,Y_,Y }, where X ± := ±(|00) ± |ll>)«00j ± (11|), 
X := I-X+-X- and Y± := |(|00)±*|ll))«00|±t<ll|) ) 
Y :=I-Y+-Y-. 

3) Then T = \[{X + - X_) + i(Y + - Y_)}; hence in 
this case, the evaluation of A requires the experimental 
specification of two independent probabilities, because 
for the given initial states the outcomes X a = Y a = 
|01> <01 1 + 1 10) <10| cannot appear. 

The following examples will demonstrate that if one 
uses the observed probabilities in a nonlinear way (i.e., by 
taking a multivariate function / = f(Xi, . . . , X n )), there 
are even simpler experiments that reveal the value of A. 
The X\ , . . . , X n are linear operators whose expectation 
values TrQ\Xj are experimentally measured. For exam- 
ple, initialize the input state to be |+) = (|0) + |l))/v2- 
The measurement of a x — |0)(1| + |1)(0| and <r y = 
— i(|0)(l| — |1)(0|) results in expectation values 

x = Tr[£ x [\+)(+\}a x ], (3.5) 
y = Tr[£ x [\+)(+\}a v }, (3.6) 
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respectively. One can evaluate the rate A by using the 
nonlinear formula 

A = v^ 2 + y 2 • (3.7) 

It is clear that this procedure is experimentally much less 
demanding than the previous alternative. 

Let us change a bit our perspective and interpret the 
same setup as an experiment in which two copies of the 
channel are used in each run of the experiment. In other 
words, let us assume that the test state is a two-qubit 
state |+)(+| <8> |+)(+| entering the channel £\®£\. The 
two-output state is then measured in a factorized mea- 
surement a x ® <J y . The general process measurement 
scheme of such two-copies type leads to a statistics 

p(F) = Tv[(n x ® Q X )F] , (3.8) 

where F are the operators (elements of the PPOVM act- 
ing on a four-qubit system) describing the individual ob- 
servations. In the considered case the outcomes are as- 
sociated with the operators 

F±,± := |++)<+ + |®|±a;,±y)<±a;,±y| ) 

where | ± x), | ± y) are the eigenvectors of a x and <j yi 
respectively. In particular, 

P(F±,±) = \[1± Re(w i)] • [1 ± Im( % )] , (3.9) 

and hence the observed probabilities are nonlinear 
(quadratic) functions of the channel's parameters Woi and 
cjio. Taking into account the fact that x = Re(woi) and 
y = im(woi), we see how to play with the observed prob- 
abilities to recover Eq. (3.7). 

The message of this example is that certain nonlin- 
ear parameters can be accessed in experiments in which 
two copies of the channel are tested simultaneously. The 
question that remains is whether the value of A can also 
be observed directly by measuring a single experimen- 
tal quantity Q only, i.e., whether there exists a single 
observable Q such that A = /(Tr[(fJ A ® £l\)Q])- 

To this end let us consider yet another experiment in 
which instead of initializing the two-qubit test system 
in the state | ++)(++ |, we send the entangled vector 
state \(f> + ) (with the definition \<p±) := (|01) ± |10>)/a/2) 
through the channel £ x ® £ x . The output state reads 

J = |[|01><01| + |10)(10| + kn| 2 (|01}<10| + |10)<01|)] , 

where (|01)(10| + |10)(01|)/2 = \<P + ){<t>+\-\<f>-){<t>-\ =: Z 
is a Hermitian operator and TrZ(|01)(01| + |10)(10|) = 0. 
That is, Z has the expectation value 

Tt[Zlj'] = ||^oi| 2 , (3.10) 

and hence the statistics of the nonlocal, two-valued 
measurement associated with the Hermitian operator Z 
uniquely and directly determines the parameter A = 



\/2(Z). This is exactly the direct method for estimat- 
ing the inverse decoherence rate that we were looking 
for. This measurement determines only a single indepen- 
dent parameter which can be immediately and uniquely 
transformed into the value of A. 

Let us note that the process POVM considered above 
is described by F± — \<fi + )((f) + \ <g> \(j)±)((j)±\ and F = 
\<t>+){<t>+\ ® (|00)(00| + 1 1 1> <1 1 1) being the output with 
vanishing probability of appearance. It is straightforward 
to verify that 

Tr[(£lA ® Q\)Q] = ^|^oi| 2 , (3.11) 

where Q = F + — F- = \<t> + ){<t>+\®Z. In conclusion, given 
that the decoherence basis is known, the value of A can 
be directly accessed and identified in a single, effectively 
two- valued, measurement. 

B. Unknown decoherence basis 

In this section we will ask the same question about the 
inverse decoherence rate but relax the assumption that 
the decoherence basis is known. Since the previous prob- 
lem is a special case of this one, it is immediately evident 
that there can be no single-copy experiment which mea- 
sures the value of A directly. Therefore, let us focus on 
the existence of direct two-copy experiments. Our goal is 
to identify the observable Q that one needs to measure 
in order to access the value of A directly. 

Let us start with a setup in which no ancilla is involved, 
i.e., the test state g € S(H passes the channel £\ ® 
£\, and we measure the expected value of a Hermitian 
operator M defined onH<E)H. Clearly, the scheme should 
be basis-independent, meaning that the operators g and 
M must both be U ® {/-invariant (have the same form in 
any basis of %). This implies that u = qg + + (1 — q)g~, 
where g± :~ P±/d±, P± are projections onto symmetric 
and antisymmetric subspaces oi'H'Si'H, respectively, and 
d± := d(d ± l)/2 are the dimensions of the subspaces 
spanned by P± . This is the well-known family of Werner 
states. Let us stress that in this section we assume d = 2, 
i.e., P+ = |00)(00| + |11)(11| + |0+X0+| and P_ = p_ = 
\(j>-){<j)-\. Similarly, any Hermitian operator M = al+bS 
(with a, b real), where S := P + — P- is the swap operator, 
is U ® U invariant. 

Let us recall once more that the expressions for P± 
have the same form in any basis of H. Therefore, 

(£ x ®£x)[P+} =A + B + X 2 Z, (3.12) 
{£ X ®£ X )[P-]=B-\ 2 Z, (3.13) 

where we used 

A:= |00)(00| + |11)(11|, (3.14) 

S := |(|01><01| + |10><10|) , (3.15) 

and the Hermitian operator Z = (|01)(10| + |10)(01|)/2 
is defined as in the previous section. Since S = A + Z = 
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|00)(00| + |11)(11| + \4>+)(<f>+\ - \4>-){<t>-l we find that 
Tv[SA] = 2, Tt[SZ] = 1 and Ti[SB] = 0. Moreover, 
TtA = 2, TtB = 1 and TrZ = 0. Using all these identities 
we obtain 

(M) = T r {(£ x ®£ x ){g}M} 

= !"&[M(j1 + B + A 2 Z)] + (1 - q)Tr[M(B - X 2 Z)} 

= | (2a + 2b + a + bX 2 ) + (1 - q)(a - b\ 2 ) 



a+-[2q + X 2 (Aq-3)]b. 



(3.16) 



Setting a = 0, b = 1 (i.e., M = S) we obtain the direct 
estimation formula for the inverse decoherence rate in its 
simplest form 



A = 



/3(S)-2g 
4a -3 



(3.17) 



Let us stress that the basis-dependent operators 
A, B, Z need not be measured; only the basis- 
independent swap operator S has to be measured. The 
whole experiment, depicted in Fig. 1, is indeed basis- 
independent and directly determines the value of the non- 
linear parameter A, as desired. 




FIG. 1. Direct inverse decoherence rate measurement, con- 
sisting of a preparation of a Werner state, two applications of 
the decoherence channel, and a SWAP measurement. 



IV. AVERAGE INVERSE DECOHERENCE 
RATE FOR QUDITS 

The case of d-dimensional quantum systems (qudits) is 
a lot more complicated because each pair of decoherence 
basis elements \j), \k) is associated with its own inverse 
decoherence rate Xjk — \wjk\- Let us perform the same 
experiment as in the qubit case (see Fig. 1). Suppose 
the initial test state is 



g = qg+ + (I - q)g- 



(4.1) 



where, as before, g± = P±/d±. The procedure reads as 
follows: Apply the channel £ x ®£ Xl where A now denotes 
all the inverse decoherence rates Xjk- And then measure 
the expectation value of the swap operator S. 

Let us define for j ^ k the following (mutually orthog- 
onal) unit vectors 



i>±,ik :=^{\jk)±\kj)) 



(4.2) 



Then 

P+ = zZ \33)UJ\ + £ \<f>+,ik)(<f>+,jk\ , (4.3) 

3 j<k 

P-=/Z^-^^-^\- ( 4 - 4 ) 

j<k 

Since = ujjk\j}(k\ and cojj = 1, we have that 

(£ x ®£ x )[P + ]=J2\jj){jj\+ (4.5) 

3 

+ l 1 £{\3k)W+^k\3k)(kj\) , 
(£x ® £x)[P-\ = W (\jk)(jk\ X%\jk){kj\) . (4.6) 



Let us note that \cJjk\ = ^jk = ^kj = \<^kj\ an d hence 

Tr (S(£ x ® £x)[P+\) =d + + J2 X %, (4-7) 

j<k 

Tr(S(£ x ®£ x )[P-]) = 0-J2x 2 jk . (4.8) 

j<k 

Now it is straightforward to see that the measured aver- 
age value of the swap operator in the considered experi- 
ment equals 

(S) =Ty(S(£ x ®£ x ){q}) 

= f c+E^-^E**- (4-9) 



j<k 



(L 



j<k 



Since cL = d(d— l)/2 coincides with the number of pairs 
jk we may define the average decoherence rate A 2 := 

rS,<t x % and obtain 

(S) = \y{qd 2 -d+) + qd]/d + ; (4.10) 

so the direct estimation formula for A 2 reads 

d+(S)-qd (d+l)(S)-2g 
qd 2 -d+ 2qd - (d + 1) ' [ ' ' 



V. DOUBLE-COMMUTATOR MASTER 
EQUATION 

In case of qubits it was shown [9] that master equations 
inducing pure decoherence channels are of the double- 
commutator form 

d g = C(g) = ~[g, H] - ^[[g, H],H] , (5.1) 



dt 



2r 7 l 



where H is a Hcrmitian operator. Let us call this oper- 
ator "Hamiltonian" and its eigenvalues "eigenenergies" . 
The channels £ = e ct are pure decoherence channels, 
and the eigenbasis of H represents the decoherence basis. 
This master equation has perfect sense in any dimension 
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(c.f. [10, 11]), and it is natural to ask whether our pre- 
vious analysis can help somehow to directly measure the 
parameter 7. 

Suppose hj are the eigenenergies of H, and denote their 
differences by Ajk := hj — hk- Then the master equation 
(5.1) reads 



c (p)= E (!-^) 

j,k=0 



Ajk 1 . Aj fc 



h V 2h l 



Pjk\j)(k\- (5-2) 



Setting ujjk — \jke lVik in Eq. (2.1), it follows from £ = e G 
that the generator of a general pure decoherence process 
has the form 

d 

9{P)= E ( 1 - S jk)( lnX ik+i<P3k)Pjk\j)(k\- (5.3) 
j,k=o 

Comparing the real and imaginary parts of Eqs. (5.2) and 
(5.3) we have 



and 



ifjk = Ajk/h 



X 3k = exp ( —i- 



2h z j 



(5.4) 



(5.5) 



With the help of Eq. (5.5) one can express the average 
inverse decoherence rate as a function of 7, 



A 2 



4 

h 2 j 



(5.6) 



The question is whether it is possible to invert this for- 
mula. If this can be done then the experimental obser- 
vation of A 2 via measurement of (S) (as described in the 
previous section) directly determines 7 = /(A 2 ). Clearly, 
if complete information about the differences Ajk of the 
eigenenergies of H is available (with the eigenbasis re- 
maining unspecified), Eq. (5.6) specifies the parameter 
7 implicitly. Unfortunately, we are not able to explic- 
itly perform the inversion, nor to specify the conditions 
under which such an inversion is possible, for an arbi- 
trary Hamiltonian. Let us note that once the parameter 
7 is specified, the knowledge of Ajk would allow us to 
determine all inverse decoherence rates Xjk- 

A special case in which the inversion is in fact possible 
is a Hamiltonian which is equally gapped (for instance, 
that of a linear harmonic oscillator) , for which 



Ajk = (j - k)A , 



(5.7) 



where A is the elementary energy gap. The absolute 
values of the differences \j— k\ take values n = 1, . . . , d— 1, 
and each of these appears exactly N n = d — n times in 
the sum present in Eq. (5.6). Therefore, Eq. (5.6) can be 
rewritten in the form 



A2 = ^E( rf --) 



K n 



(5.8) 



where K := exp [-A 2 /(^ 2 7)] . 
formula for gamma reads 



where 



K r := Root r 



h 2 lnK r 



d-l 



The direct estimation 



(5.9) 



x 2 d- -J2( d ~ n ) Rr ' 



n=l 



(5.10) 



is the rth root of the above polynomial in the variable 
K with r = 1 for even Hilbert space dimension d and 
r = 2 for odd d. The ordering of roots is such that real 
roots precede complex ones, with the real roots ordered 
by size. We have verified numerically for all d up to 32, 
that the rth root, Eq. (5.10), is the unique positive real 
solution of Eq. (5.8). 

Morevoer, the identity K r — exp [— A 2 /(7i 2 7)] and 
Eqs. (5.5) and (5.7) enable us to write down the esti- 
mation formula for the individual off-diagonal elements 



(5.11) 



n=l 



Let us note that in this case the knowledge of the actual 
value of the gap A is not necessary. 

In summary, the method of direct observation of the 
average inverse decoherence rate A 2 can be employed [cf. 
Eq. (5.6)] to estimate the decoherence decay parameter 7 
which occurs in the double-commutator master equation 
given in Eq. (5.1), provided that the additional knowl- 
edge on the structure of energy gaps (spectrum) is avail- 
able. This information can then also be used to determine 
directly all individual inverse decoherence rates Xjk ■ 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we proposed an experimental setup that 
can directly determine the value of the average inverse 
decoherence rate, provided that the channel under con- 
sideration is a so-called pure decoherence channel. The 
experiment consists of the following steps: 

1. Prepare a bipartite Werner state with q 7^ d + /d 2 
[see Eq.(4.1)]. 

2. Apply the channel £\ (Si £\- 

3. Perform the SWAP measurement and record the 
expectation value. 

4. Apply the estimation formula (4.11). 

Let us recall that the SWAP measurement is a two- valued 
measurement, and hence by definition, it can fix only 
one parameter of £\. Therefore, no other information 
about decoherence beyond its average rate is revealed. 
For instance, the decoherence basis remains as unknown 
as it was at the beginning. 
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Let us stress that no complicated system source is 
needed in order to accomplish this measurement. In fact, 
any state g can be turned into a Werner state by a so- 
called twirling transformation 



u q = J dUU ® UqU* ® . 



(6.1) 



Moreover, instead of a source of bipartite systems we can 
employ two copies rj®rj produced by the same single par- 
ticle source. Applying the same random unitary channel 
on each of the copies will result in a Werner state with 
some initial value of q (see Fig. 2). 




FIG. 2. Preparation of a two-partite Werner state using a 
single-particle source and sampling of a random unitary chan- 
nel 



This value is related to the purity of the single particle 
source by 



q= ^(1 + TV). 



(6.2) 



This is the only information about the test state that 
we really need, and it is well known how to measure the 
purity directly [20, 21]: Take two copies, perform the 
SWAP measurement and relate Trr/ 2 = Tr Sij (8177= (S). 
If one uses this way of preparing the Werner state, the 
accessible values of q are limited by d + /d 2 < q < 1. 
Furthermore, let us note that the mixing procedure pre- 
supposes that no information on the actual sequence of 
the random Us is used in the estimation. Otherwise, the 
procedure could not be seen as a preparation of a Werner 
state. 

The key operation in the preparation procedure is the 
twirling. It is experimentally difficult to really sample 
from the Haar measure over unitary channels. However, 



this is not really needed, as the integration can be re- 
placed by an average over only a discrete number of uni- 
tary channels. This construction is known as a unitary 
fc-design [17-19]. In our particular case we are inter- 
ested in a so-called 2-design, which is any probability 
distribution {pi, . . . ,p n } distributed over unitary opera- 
tors {Ui, . . . , U n } such that for all operators A, 



dU(U® 2 )X(U®y =Y,Pj(Uj ® Ui)X{U] ® U}) . 
jeJ 



In particular [17], for a d-dimensional quantum system 
the uniform 2-design (i.e. pj = 1/n) consists of n > 
d 4 -2d 2 + 2 unitaries. 

Let us summarize which resources are sufficient for di- 
rect observation of the inverse decoherence rate: The ex- 
perimenter needs a single particle source of (not com- 
pletely mixed) states; he must be able to implement a 
finite number of (undisclosed) unitaries forming the sup- 
port of a 2-design; and finally, he must be able to perform 
the SWAP measurement. 

In conclusion, we have addressed the problem of direct 
experimental measurement of a specific nonlinear chan- 
nel parameter, the average inverse decoherence rate. We 
have shown that experiments which collect statistics on 
two copies of the channel perfectly serve this purpose. 
Since many interesting channel properties (quantum ca- 
pacities, rates, etc.) are not linear and cannot be accessed 
directly in single-copy experiments, it is of practical im- 
portance to understand more deeply the theory of direct 
estimation of nonlinear channel parameters. The present 
paper represents an important case study in this regard. 
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